Towards a wave-extraction method for numerical relativity. 
IV. Testing the quasi-Kinnersley method in the Bondi— Sachs framework 
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We present a numerical study of the evolution of a non-linearly disturbed blaclt hole described 
by the Bondi-Sachs metric, for which the outgoing gravitational waves can readily be found using 
the news function. We compare the gravitational wave output obtained with the use of the news 
function in the Bondi-Sachs framework, with that obtained from the Weyl scalars, where the latter 
are evaluated in a quasi-Kinnersley tetrad. The latter method has the advantage of being appli- 
cable to any formulation of Einstein's equations — including the ADM formulation and its various 
descendants — in addition to being robust. Using the non-linearly disturbed Bondi-Sachs black hole 
as a test-bed, we show that the two approaches give wave-extraction results which are in very good 
agreement. When wave extraction through the Weyl scalars is done in a non quasi-Kinnersley tetrad, 
the results are markedly different from those obtained using the news function. 

PACS numbers: 04.25.Dm, 04.30.Db, 04.70.Bw, 95.30.Sf, 97.60.Lf 
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I. INTRODUCTION 

Gravitational-wave detection has been gaining much 
interest over the last decades. Much effort has been cor- 
respondingly made in modelling possible sources. Inter- 
esting examples of sources of gravitational waves are, e.g. , 
binary systems of merging black holes, spiralling systems 
of two neutron stars or coalescing black hole-neutron star 
binaries. Analytical tools able to investigate the dynam- 
ics of such sources when the merging takes place do not 
exist, due to the strong non-linearity of the problem. Nu- 
merical simulations are therefore invaluable in order to 
extract information about the gravitational-wave signal 
emitted. Typically, numerical relativity studies are done 
in three stages; First, one specifies initial data that cor- 
respond to the physical system of interest, and that sat- 
isfy certain constraint equations. Next, one evolves these 
initial data numerically, using the evolution equations 
(with or without enforcing the constraints), and finally, 
one needs to interpret the results of the simulation and 
extract the relevant physics thereof. This paper — like 
its prequels P, 0, 13 — is concerned with this last stage, 
namely, with the problem of wave extraction. 

In order to interpret the results of numerical relativ- 
ity simulations — at least within the context of a specific 
formulation of Einstein's equations — a useful approach 
is that based on the characteristic initial value problem, 
originally introduced by Bondi and Sachs 0, Q. The 
characteristic initial value problem has been extensively 
used in numerical relativity, for spherically symmetric 
systems lEHJl H EI IH IM, axisymmctric sys- 
tems [lllESlS and 3D systems^J, Ilkll9. .2QJ. This 
article is based on the results found in [13 — where the 
characteristic initial value problem, using Bondi coordi- 



nates, is used in axisymmetry to study a non-linearly 
disturbed non-rotating black hole metric. The non-linear 
response of the Schwarzschild black hole to the gravita- 
tional perturbation is embodied in a superposition of an- 
gular harmonics propagating outside the source. All the 
information concerning the angular harmonics and the 
energy radiated are easily derived by the evolved quan- 
tities. In fact, in this particular case one can identify a 
news function i.e., a function which embodies the in- 
formation about the gravitational-radiation energy emit- 
ted m mm mil mm. 

It is not currently possible to have at hand a quantity 
such as the news function when using other numerical 
approaches like, specifically, the widely used 3-1-1 de- 
composition of Einstein's equations [23 • Indeed, one 
of the outstanding problems of numerical relativity is 
that of wave extraction, i.e., the problem of how to ex- 
tract the outgoing gravitational waves from the results 
of numerical simulations. In currently available meth- 
ods, various approximations are applied to determine 
the gravitational-wave emission of isolated sources. One 
of the simplest approaches applies the quadrupole for- 
mula (that strictly speaking is valid for weak gravita- 
tional fields and slow motions) mi • This approach 
has been used effectively, e.g., in models of stellar col- 
lapse [33 ■ More sophisticated approaches use the Mon- 
crief formalism [sills^ to extract first-order gauge invari- 
ant variables from a spacetime which is assumed to be a 
perturbation of a Schwarzschild background at large dis- 
tances js^, HE mi ■ The strength of this approach is 
that it is gauge invariant, i.e., the information extracted 
is related to the physics of the system, and not to the 
coordinates used. Specifically, it avoids the pitfall of 
misidentifying gauge degrees of freedom as gravitational 
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waves. [E.g., in the Lorenz gauge all degrees of freedom, 
including the (residual-) gauge ones, travel at the same 
speed — the speed of light — such that they might be con- 
fused with physical waves.] This procedure is usually per- 
formed under the assumption that the underlying gauge, 
i.e., the particular choice of the spacetime coordinate sys- 
tem, leads to a metric which is asymptotically Minkowski 
in its standard form, which is indeed the case for the most 
commonly used gauges in simulations of isolated systems. 
The gravitational waveform is determined by integrating 
metric components over a coordinate sphere at some ap- 
propriately large distance from the central source, and 
then subtracting the spherical part of the field (which is 
non-radiative). However, as we have already mentioned, 
such techniques are well defined when the background 
metric is assumed to be Schwarzschild, while their appli- 
cation to the more generic Kerr background metric can 
at best be intended as a very crude approximation. In ad- 
dition, in the typical numerical relativity simulation, one 
does not usually have information about the mass and 
spin angular momentum of the eventual quiescent black 
hole, and no prescription is currently known to uniquely 
separate a background from a perturbation. Moreover, 
once the elimination of the spherical background is per- 
formed, what is left does not necessarily satisfy the per- 
turbative field equations, and may, in fact, be quite large 

M 

Another important approach is that of Cauchy- 
charactestic matching (CCM) schemes The pecu- 

liarity of such schemes is that they use a characteristic 
Bondi-Sachs approach to study the numerical space-time 
far from the source, where the fields are weak and the 
probability to form caustics, which would make the code 
crash, is limited. In this way, using the notion of the 
Bondi news function, it is possible to extract easily the 
gravitational wave information. In the strong field part 
of the computational domain, instead, a usual Cauchy 
foliation is used, so that the problem of caustic forma- 
tion is irrelevant. The CCM schemes have been used 
successfully to simulate cylindrically symmetric vacuum 
space-times ^3 or to study the Einstein-Klein-Gordon 
system with spherical symmetry ^ fully 3D appli- 
cation of the CCM scheme is, however, still unavailable. 

A final approach worth mentioning here aimed at wave 
extraction is the one involving the Bel-Robinson vector 
|38j | , which can be considered a generalization to general 
relativity of the Poynting vector defined in electromag- 
netism. However its connection with the radiative de- 
grees of freedom is still not entirely clear. 

A novel approach has been suggested recently pi 13.13^: 
one extracts information about the gravitational radia- 
tion through quantities that are gauge and background- 
independent. One such quantity is the Beetle-Bur ko 
scalar, which is also tetrad independent. Specifically, 
no matter how one chooses to separate the perturba- 
tion from the background, the Beetle-Burko scalar re- 
mains unch ang ed. However, as pointed out in various 
contexts [3, llfl, Ell , the physical meaning of the Beetle- 



Burko scalar is non trivial. For example, in the stationary 
spacetime of a rotating neutron star its non-zero value 
is due to the deviation of the quadrupole from that of 
Kerr while clearly no radiation is present. Thus 
the Beetle-Burko scalar awaits further study. At any 
case, the Beetle-Burko scalar — while including informa- 
tion only on the radiative degrees of freedom (when the 
notion of radiation is defined unambiguously) — describes 
the latter only partially. To obtain a full description of 
the radiative degrees of freedom, it is therefore desirable 
to consider an approach in which one calculates quan- 
tities, whose physical interpretation is more straightfor- 
ward; however, the price to pay for these advantages is 
that it is harder to obtain such quantities uniquely: in 
fact, one still needs to break the spin/boost symmetry in 
a useful way. 

This approach, which is the basis of this paper, is that 
of using the Weyl scalars — which, under certain assump- 
tions, isolate the radiative degrees of freedom from the 
background and gauge ones [a^E3jE3 — for wave extrac- 
tion. Teukolsky %4\ showed that, choosing a particular 
tetrad, namely the Kinnersley tetrad |45| . to calculate 
the Weyl scalars, it is possible to associate the Weyl 
scalar ^4 with the outgoing gravitational radiation for 
a perturbed Kerr space-time. Other authors have sug- 
gested to use the Weyl scalars for wave extraction, most 
recently in (46| . or to explore their relation with met- 
ric perturbations However, the Weyl scalars de- 
pend on the choice of tetrad. Specifically, performing 
null rotations on the basis vectors of the null tetrad one 
can change the values of the Weyl scalars. (Recall, that 
one of the advantages of the Beetle-Burko scalar is that 
it is tetrad independent.) In fact, extracting the Weyl 
scalar ^4 is meaningless, unless one also describes how 
to construct the tetrad to which it corresponds. In most 
tetrads, the Weyl scalar ^'4 mixes the information of 
the outgoing radiation with other information, includ- 
ing gauge degrees of freedom. In our proposal the use 
of the Weyl scalars is intimately related to the construc- 
tion of the tetrad in which they are to be calculated — the 
quasi-Kinnersley tetrad — and therein lies its strength. 

With the aim of using '^4 to extract information from 
simulations about the gravitational radiation output, 
starting with a transverse condition "ifi = '^3 = 0, recent 
work Q, H, 111 has addressed the problem of computing 
Weyl scalars in a tetrad which will eventually converge 
to the Kinnersley tetrad. This tetrad has been dubbed 
the quasi-Kinnersley tetrad, and work is still in progress 
in order to uniquely identify it for a general metric. Up to 
now, it is possible, for a general spacetime, to identify a 
class of tetrads, namely the quasi-Kinnersley /rame 
with the property that in this frame the radiative degrees 
of freedom (when and where the notion of radiation is 
unambiguous) are completely separated from the back- 
ground ones. However, work is still in progress to identify 
the quasi-Kinnersley tetrad out of this frame. The dif- 
ficuly in doing so is related to the following property of 
the tetrad members of the quasi-Kinnersley frame: they 
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are all connected through type III ("spin/boost") rota- 
tions, and the spin/boost symmetry needs to be broken 
before the Weyl scalar ^4 can be extracted in the quasi- 
Kinnersley tetrad. The quasi-Kinnersley frame is there- 
fore a two-parameter family of tetrads, and the value of 
^1^4 (and that of 5*0) depends on the choice of the tetrad 
member of the frame. Notably, the Beetle-Burko scalar is 
invariant under type III rotations. That is, all the tetrad 
members of the quasi-Kinnersley frame share the same 
Beetle-Burko scalar. We avoid the spin/boost symmetry 
breaking difficulty in the present paper by applying an ad 
hoc technique to find the quasi-Kinnersley tetrad. This 
ad hoc technique allows us to obtain the Weyl scalar '^4 
in the quasi-Kinnersley tetrad, where its interpretation 
according to the gravitational compass is readily avail- 
able. 

Specifically, in this paper we use the Bondi-Sachs for- 
malism i, as it turns out that in this special case we can 
identify a (non-transverse) quasi-Kinnersley tetrad in a 
simple way, and to compute the Weyl scalars directly. 
The aim of this work is thus to demonstrate — in the con- 
text of this practical numerical example — the applicabil- 
ity, and necessity, of the quasi-Kinnersley tetrad method 
in using the Weyl scalars as wave extraction tools. 

The article is organized as follows: Section |nl intro- 
duces our physical scenario, and Section lllll describes our 
Weyl scalars computation. In Section Hvl we present the 
expected result which links the Bondi news function to 
^1^4. Finally results and conclusions are presented in Sec- 
tion E| 



II. THE BONDI PROBLEM 

The numerical scenario we are studying is that of a 
non-linearly perturbed^ Schwarzschild black hole using 
an ingoing null-cone foliation of the space-time. 

We set up our system of coordinates as follows: a time- 
like geodesic is the origin of our coordinate system, pho- 
tons are travelling from the origin in all directions, their 
trajectories forming null hypcrsurfaces. The hypersur- 
face foliation is labelled by the coordinate v. As r co- 
ordinate we choose a luminosity distance, such that the 
two-surfaces of constant r and u have area 47rr^ . Finally, 
each null geodesic in the hypersurface is labelled by the 
two angular variables 9 and (f). We will restrict our at- 
tention to an axisymmetric space-time such that ^ is a 
Killing vector. Having chosen these variables, the Bondi 
metric in ingoing coordinates reads 



^ We want to point out here that the expression "perturbed" could 
be misleading in this context, as it might suggest we are assum- 
ing some kind of approximation. Our numerical simulations are 
instead fully non-linear evolutions of Einstein's equations in the 
Bondi-Sachs formulation. 



M 



dv' 



2e^'^dvdr - 2Ur^e^^dvde 



(1) 



where M, U, (3,j are unknown functions of the coordi- 
nates {v, r, 9). Within this framework, the Einstein equa- 
tions decompose into three hypersurface equations and 
one evolution equation, as given below in symbolic nota- 
tion 



□ (2)^ = n^{M,(3,U,-f), 



(2a) 
(2b) 
(2c) 
(2d) 



where d'^^) is a 2-dimensional wave operator, ip = rj, 
and the various H symbols are functions of the Bondi 
variables. We will write here only the expression for Hp, 
the simplest of all the functions, which is given by 



n 



/3 



-r7 



(3) 



We will use this expression to describe the numerical al- 
gorithm. For exhaustive description of the system ^ we 
refer to [13 and 

The structure of Eq. (01 establishes a natural hierar- 
chy in integrating them. By setting the initial value for 
the function 7 on the initial hypersurface (in addition to 
four free parameters), it is possible to integrate Eq. (|2b|) 
to obtain f3, then, having both /3 and 7, Eq. flc^ can 
be integrated to obtain U and finally M can be derived 
by integrating Eq. Ij2d|l . At this point we have all the 
metric functions on the initial hypersurface and we can 
integrate Eq. H2a|) to obtain 7 on the next hypersurface, 
and the procedure is iterated. The integration of Eq. (O 
introduces constants of integration, which we set to zero 
(Bondi frame 14^1. which corresponds to having an 
asymptotically inertial frame. 

The metric introduced in Eq. ^ describes a static 
Schwarzschild black hole if we set, in the Bondi frame, 
all the functions except M to zero everywhere in the do- 
main. M is chosen to be the Schwarzschild mass Mq of 
the black hole. Besides, outgoing gravitational radiation 
as a perturbation is introduced by the function 7; it turns 
out that 7 is a spin-2 field and is actually related to the 
radiative degree of freedom. Choosing an initial shape 
for 7 means in practice choosing the initial profile of out- 
going gravitational waves. More specifically, the initial 
data are chosen in the following way: 

• We recover the background Schwarzschild geome- 
try by setting -f — (3 — U = over the whole 
computational domain, while we set M = Mq 
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FIG. 1: Diagram of the Bondi algorithm: the foliation is built up using ingoing light cones that emanate from an external 
worldtube W. We fix the geometry at the world tube to be that of a single Schwarzschild black hole. We non-linearly perturb 
the space-time by setting a non- vanishing value for the Bondi function 7. Such disturbance propagates outwards in the radiative 
region, after scattering ofi^ the black hole. 



• We set up the initial data for a gravitational wave 
outgoing pulse by choosing a gaussian shape (with 
parameters Tc and cr) for the function 7 (the choice 
of the gaussian shape is dictated by requiring the 
function to vanish at the outer boundary of the 
grid): 



7(r,0) = -=_e-'^y2,„(0), (4) 
V27rcr 

where A is the amplitude of the perturbation, and 
Yim is the spherical harmonic of spin 2. In our nu- 
merical simulations we will set two types of initial 
data: the first one with I = 2 and m = to get 
a pure quadrupole outgoing perturbation, while in 
the second one we will set / = 3 and m = 0. 

The integration of the hypersurface equations leads to 
the problem of the gauge freedom in choosing the integra- 
tion constants. This apparent freedom is fixed by choos- 
ing outer boundary conditions on our numerical grid. As 
depicted in Fig. we in fact fix the metric at the outer 
boundary, i.e. on the worldtube W, to be that of a 
Schwarzschild black hole. This automatically fixes the 
integration constants to be for 7, U and /3, and Mq for 
M (Bondi frame). We point out here that such a bound- 
ary condition is well posed only for simulations which are 
limited in time, so that no relevant outgoing gravitational 
flow has crossed the worldtube. More information about 
the evolution routine can be found in |,l&i49t .i50l|. 



III. WEYL SCALARS 

Once we have the numerically computed metric for the 
evolved space-time we can derive the Newman-Penrose 
quantities we need. The Weyl scalars are 



*0 = 




(5a) 


*1 = 




(5b) 


*2 = 


^Cabcdrm''m''n^ 


(5c) 


*3 = 


~Cabcdl''n''m''n\ 


(5d) 


^1^4 = 


-Cabcdn'^ml'rffh'^, 


(5e) 



where n", m" and m° are the Newman-Penrose null 
vectors. The five scalars defined in Eq. are of course 
coordinate independent, but they do depend on the par- 
ticular tetrad choice. 

We calculate the scalars in a quasi-Kinnersley tetrad, 
i.e., in a tetrad that converges to the Kinnersley tetrad 
when space-time settles down to that of an unperturbed 
black hole. In P, a procedure to find the quasi- 
Kinnersley frame in a background independent way is 
given by looking at transverse frames, i.e., those frames 
where 4'i and '^^ vanish. In the following we use the word 
frame to indicate an equivalence class of tetrads which 
are connected by a type III spin/boost tetrad transfor- 
mation. Fixing the right quasi-Kinnersley tetrad means 
choosing the tetrad in the quasi-Kinnersley frame which 
shows the right radial behavior for and \E'4 according 
to the peeling-off theorem. That is, we seek the tetrad in 
which these two Weyl scalars peel off correctly. A back- 
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ground independent procedure to single out the tetrad 
out of the frame currently needs further investigation. 

In the Bondi-Sachs framework the identification of 
a quasi-Kinnersley tetrad is simple, and does not need 
to use the notion of transverse frames. The main rea- 
son for this proprety is due to the asymptotic knowl- 
edge of the Bondi functions when spacetime approaches 
Schwarzschild: in fact 7, U and /3 tend to zero, while M 
tends to the Schwarzschild mass Mq of the black hole (for 
further details see the next section). 

This situation is much different from the typical sit- 
uation in numerical relativity simulations, for which 
the wave-extraction methods needs to be background- 
independent. If we assume to be in the Schwarzschild 
limit, the background Kinnersley tetrad chosen by 
Teukolsky ,44] in the perturbative scenario would look, 
using our coordinates {v,r,6,(j)), as 



2r 



r-2M 



,1,0,0 



„ r ~ 2M 
0, ,0,0 



2r 



0,0, 



V2r' V^rsm0 



(6a) 
(6b) 
(6c) 



This tetrad has been chosen by letting the and 
vector coincide with the repeated principal null directions 
of the Schwarzschild space-time. Such a condition fixes 
a frame, i.e. a set of tetrads connected by a type III 
rotation. The type III rotation parameter is then fixed 
by setting the spin coefficient e to be vanishing. Eq. © 
can be used to find the general expression for the tetrad 
in the full Bondi formalism. Using the asymptotic values 
of the Bondi functions, we can write down the expression 
for a general tetrad for the Bondi metric, whose vectors 
converge to the null vectors written in Eq. ^ in the 
Schwarzschild limit. The result is given by 



[(1 - 2M/r) e4/3 - {72^2g2(7+/3)- 



0,0 



[(1 - 2M/r) e2/3 - UVe^^] 



0, 



r[/e(7-2/3) 1 



V2 ' V2reT' %/2r sin 6*6-7 



(7a) 
(7b) 

(7c) 



Henceforth we will denote this first tetrad choice, which 
is supposed to be the successful one, as the tetrad Ti. 
It is worth pointing out that the tetrad Ti is not trans- 
verse, i.e. ^1 and ^"3 are not vanishing; Nevertheless it 
satisfies the requirements needed of a quasi-Kinnersley 
tetrad. (The quasi-Kinnersley tetrad does not have to 
be transverse, although it does need to be asymptotically 
transverse.) 



As explained above, the choice of Ti has been driven by 
the form of the Kinnersley tetrad expressed in Eq. for 
the unperturbed black hole. However, this is not always 
possible in general where tetrad choices are dictated by 
different criteria. A straightforward example in this case 
would be that of basing the tetrad on null vectors directly 
derived from the metric. For instance, in the specific 
example of the metric an apparently natural choice 
of the tetrad, obtained with the algebraic manipulation 
packages Maple and GRTensor, is the following: 



[0,-6-4^0,0] 



[{l-2M/r)e'f'-UW] 

^ 1 n 1^1^ 



0, 



y/2 ' \/2reT ' V2r sin 9e 



(8a) 

(8b) 
(8c) 



The reason why this tetrad choice looks more natural 
than the first one is related to the fact that packages 
like GRTensor construct this tetrad starting from the 
vector, which is assumed to be lying on the null foliation, 
leading to the expression £^ = (5^o- 



The contravariant 



components are then given by Eq. H8a|l . Once is fixed, 
the other tetrad vector expressions are found by imposing 
the normalization conditions between the vectors in the 
Ncwman-Pcnrose formalism. 

The tetrad in Eq. ||SJ will be hereafter referred to as 
tetrad 72. We will show that this different tetrad choice 
leads to results which, although equivalent from a quali- 
tative point of view, are different than those obtained in 
the quasi-Kinnersley tetrad Ti . 

In the numerical results that we are going to present in 
the next sections we have calculated the Weyl scalars in 
the two tetrads presented above, to show the reliability of 
quasi-Kinnersly tetrad method and an example, through 
a complete comparison, of the different results that we 
might obtain in a numerical simulation in doing wave 
extraction using Weyl scalars in different tetrads. 



IV. THE LINEAR REGIME 

In the linear regime both the Bondi and the Newman- 
Penrose formalisms define quantities which provide infor- 
mation about gravitational waves. We will discuss here 
briefiy such definitions in order to get a correspondence 
between the two approaches, which will be tested numer- 
ically in the following section. 

In the Bondi formalism the initial assumption is that 
the space-time is asymtotically flat, which leads to the 
following expansion for the function 7 at null infinity: 



(9) 
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Here we assume to be in the Bondi frame i.e. we set 
the integration constant K to be zero. By integrating the 
hypersurface equations for the other Bondi functions we 
can get their radial expansion at nuU infinity. Such inte- 
grations in general introduce other integration constants 
but, again, we assume that in our frame those constants 
are vanishing, ending up with the following expressions 
for the remaining Bondi functions f3 | : 



/3 

U 
M 



c 



Oir- 



csm 



r-^ sin^ b 
Mo + O (r-i) . 



O r- 



(10a) 

(10b) 
(10c) 



It is trivial to verify that, with this choice of integration 
constants, the space-time is asymptotically flat. It is pos- 
sible to define at null infinity a notion of energy, which 
leads to the result found by Bondi 



E ^ — (/) A/ sin 

in 



(11) 



and in addition the energy fiux per unit solid angle, which 
is given by 



The components of the Riemann tensor in Eq. (|14() can 
then be related to the transverse-traceless gauge terms, 

using Rvav/3 — — \ Q^^'^ 1 which leads to the result 



(*4) 



qKT 



d^hee 



I - 



(15) 



This relation between \l/4 in the quasi- Kinnersley tetrad 
and the transverse-traceless (TT) components of the per- 
turbed metric leads us to a definition of the energy emit- 
ted by simply calculating the expression of the energy 
tensor for the gravitational wave defined by 



rpGW _ ^ 

327r 



dAh 



TT\ 



(16) 



The total energy fiux is then given by the formula, 
which is assumed to hold at null infinity: 



dvdO. 



r 

T6n 



dh 



TT 



dv 



dh 



TT' 



dv 



(17) 

and, by substituting our expression in terms of ^'4 we get 
the result 



d^E 
dvdCl 



4:71 



(12) 



It is clear that the information about the energy carried 
by gravitational waves is contained in the Bondi news 
function c^v « i^J,v-, where the approximation is assumed 
to hold at large distances in the linear regime. Expressing 
Eq. H12|) in terms of 7 gives 



d^E 
dvdQ 



4:71 



(13) 



An analogous derivation can be achieved within the 
Newman-Penrose formalism. The key point is that 5*4 
can be expressed, when computed in the quasi- Kinnersley 
tetrad, directly as a function of Riemann tensor compo- 
nents, i.e. 0,13 



(*4) 



qKT 



R 



vOve 



(14) 



The hatted symbols in Eq. H14() are indicating that the 
Riemann tensor components are contracted over a tetrad 
of vectors oriented along the coordinates. In the linear 
regime however, those vectors can be assumed to be the 
basis coordinate vectors, as the perturbation is already 
contained in the Riemann tensor. For this reason we will 
omit the hatted symbols from now on, and always talk 
about coordinate components. 



d^E 
dvd^ 



r 

47T 



(*4) 



qKT 



dv 



(18) 



In our specific case of axisymmetry we expect only one 
polarization state to be present, which is made evident 
by the presence of a single news function. Correspond- 
ingly, we expect (^'4)^^^^ to have only its real part non 
vanishing. Combining Eq. H18I) with Eq. (|13|l and con- 
sidering the presence of only one polarization state, we 
finally obtain that in the linearized regime the relation 



(*4) 



qKT 



dv^ 



(19) 



must hold. This is the relation we want to verify numeri- 
cally. The minus sign comes from the negative sign given 
in Eq. H15|) . We want to stress again the attention to 
the fact that such relation is strictly true at null infinity, 
however, we expect it to be well satisfied provided we are 
at sufficiently large distances from the black hole. As it 
will be clear in the next section, this assumption turns 
out to be very well motivated. 



V. NUMERICAL RESULTS 

In this section we present numerical results for a stan- 
dard simulation. We have written a code that solves the 
Bondi equations and calculates the Weyl scalars in the 
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FIG. 2: Convergence test for the R22 component of the Ricci 
tensor. The top panel shows this component for two different 
resolutions for a radial slice on the equatorial plane. As ex- 
pected the value is converging to zero. In the bottom panel 
we have tested the second order power law of convergence by 
multiplying the 1200x120 output by a factor of four. The 
two curves now overlap perfectly, thus proving second order 
convergence. 




10 20 30 40 50 60 




-- 600x60 
— 1200x120 
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FIG. 3: (a) Top panel: the value for 'I'4 for two different 
resolutions at time v — 80. (b) Bottom panel: the value for 
\['o for two different resolutions at « = 80. Both values are 
calculated on the equatorial plane. 



two tetrads Ti and T2 . Our code makes use of the Cactus 
infrastructure |5lj |. 

We set up an initial Schwarzschild black hole, and con- 
struct an initial quadrupole perturbation on 7 using the 
expression indicated in Eq. Q). The values chosen in 
this case are A = 0.1, rg = 3 and cr = 1, although vari- 
ous tries have been performed varying these parameters, 
all leading to the same physical results. We emphasize 
that A represents the amplitude of the initial perturba- 
tion: the chosen value is such that the perturbation is 
somehow realistically small, yet large enough for the full 
non-linearity of the problem to appear clearly through 
the harmonic coupling, as we are going to show (see Figg. 
llOland im cf. also 16]). All the results presented here are 
obtained using two different resolutions, the coarser one 
having 600 points in the radial dimension and 60 points 
in the angular direction, the finer one having those val- 
ues doubled. The results which are not convergence test 
results are all obtained using the finer resolution of 1200 
points in the radial direction and 120 points in the angu- 
lar direction. 

We will first present some tests in order to verify the 
robustness of our algorithm, and then we will proceed to 
a full comparison of our results in the two approaches 
presented here. The first two following subsections will 
deal with the calculation of the Weyl scalars in the tetrad 
Ti defined in Eq. Q; we don't expect the second tetrad 
to give different results for what concerns radial fall-offs 
and convergence. Section IV CI will instead deal with the 
relation of ^'4 with the news function and, within this 
context, it is very important to show a comparison of 
results in different tetrads, to have an evident demon- 
stration of how important the choice of the right tetrad 
is, i.e. the quasi-Kinnersley tetrad, in the process of eval- 
uating the outgoing gravitational wave contribution. 



A. Convergence 

The first thing we want to test in our code is of course 
convergence. In order to do so, once the numerical vari- 
ables are computed, we have calculated independently 
the values of the Ricci tensor components, which should 
vanish in vacuum. Such components are suitable for do- 
ing convergence tests. In Fig. (|2Jl we show the value of 
the Ricci component R22 for two resolutions, the picture 
shows a radial slice of our space-time on the equatorial 
plane, for the time value v = 80. The first figure simply 
superimposes the two values obtained for the two differ- 
ent resolutions, while in the second picture we have first 
multiplied the values for the finer resolution by a factor 
of four, as expected in a second order convergence code. 

We have found similar results for the other compo- 
nents, which ensured us of the convergence of our algo- 
rithm. 
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B. Radial fall-offs of the scalars 

Figg. (PJl) and (jjja) show the numerical output for '^q 
and \l/4 for two different numerical resolutions. The out- 
puts show satisfactory convergence for ^4, but not for 
'i'o. This is because the asymptotic radial behavior for 
^0 should be r~^, as expected from the peeling-off con- 
jecture (see e.g. |^) and the linear perturbation analysis 
p4l |. and this gets completely embedded in the numeri- 
cal error. We believe that this could constitute a serious 
numerical problem in situations where the initial tetrad 
chosen for the scalars computation is not the right one, 
and a tetrad rotation is needed. The numerical error 
found in ^'o would then propagate when other quanti- 
ties, like the curvature invariants / and J, are computed, 
thus leading to meaningless results. Recall, however, that 
the curvature invariants /, J, in addition to the Coulomb 
scalar x and the Beetle-Burko scalar ^, can be found in- 
variantly and in a background-independent way which is 
also tetrad-independent, i.e., it does not require finding 
first the Weyl scalars to find I, J P, Q • 

Fig. 10^) and Qd) emphasize the radial dependence of 
^'2 and ^'4, which is highlighted very well in our simu- 
lations. The two figures show that at late times 5*2 gets 
the background contribution with the superposition of a 
wave whose radial behaviour is r~^. We have tested the 
convergence of such a wave to prove its physical mean- 
ing; this is itself a quite interesting result as we don't 
have a perturbation equation for 5*2, and it is entirely 
due to the full non-linear treatment of the problem. Of 
course, given its rapid fall-off, the wave contribution from 
^2 is negligible. ^'4 shows instead the well expected r^^ 
behaviour. 




10 20 30 40 




10 20 30 40 



r (BH units) 

FIG. 4: (a) Top panel: the value for fl2 = at two 

different times vi = 40 and V2 = 80. (b) Bottom panel: the 
value for fl4, = r>I'4 for the same couple of times vi — 40 and 
V2 = 80. 

In order to understand what is happening, we need to 
analyze further the tetrad T2, and in particular its limit 
when the space-time approaches a type D one. Using our 
well known asymptotic limits for the Bondi functions, it 
is easy to show that tetrad T2 converges, in the type D 
limit, to the tetrad 



C. Relation of ^£'4 with the Bondi news 

In this section we want to show the comparison of '^4 
with the second time derivative of 7, where 5*4 will be 
calculated in the two tetrads shown in Eq. Q and l|HJl. 
We first start with the tetrad 71 : Figg. lO and © verify 
numerically the equivalence expressed by Eq. (|19|) : it is 
clear that the two functions ^'4 and — 7,«t, are different in 
the non-linear regime but converge in the linear regime. 
In particular Fig. © shows in logarithmic scale the ab- 
solute value of their difference at time v = 80, well in the 
linear regime. This numerical result proves the generic 
assumption that '54 is related to the outgoing gravita- 
tional radiation contribution. 

As a counterexample, we show the same result when 
^'4 is computed in tetrad 72, which would actually have 
been our simplest choice hadn't we applied the concept of 
a quasi- Kinnesley tetrad. The results for this calculation 
are shown in Fig. {T)). It is evident that \I'4 does not 
get any contribution from the background, meaning that 
the tetrad we have chosen is part of the quasi- Kinnersley 
frame, however, it is evident that the result is rather 
different from that coming from the Bondi function 7. 



= [0,-1,0,0], 
r-2M 

1, 



0,0, 



2r 
1 



,0,0 



\/2r ' \/2r sin 



(20a) 
(20b) 

(20c) 



which is different from the Kinnersly tetrad used by 
Teukolsky, Eq. A simple analysis of the differences 
let us conclude that the original tetrad Eq. ijHJl can be 
obtained by first of all exchanging the two real null vec- 
tors H. and n, and then using a boost transformation of 
the type 



n, 



(21a) 
(21b) 



where A is a real parameter. It is easy to show that 
choosing A = -pz^ g^t that the new real null vectors 
coincide with the Kinnersley tetrad defined in Eq. 
It is now straightforward to understand how these dif- 
ferences affect the values of the Weyl scalars. First of 
all, exchanging € and n corresponds to exchanging ^'o 
and ^'4; this means that if we use the tetrad 72 we will 
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FIG. 5: From top to bottom: the comparison of *I'4 (dashed 
hne) calculated in tetrad 71 and —^,vv (solid line) for the 
values of va = 10, v\ = 40, i;2 = 60 and «3 = 80. 




FIG. 6; The function A = 1*4 +7,i,ul at wa 



find the outgoing radiative contribution in ^o- This com- 
pletely clarifies the result found in Fig. (jT)): it turns out 
that in this particular tetrad ^'4 is supposed to have a 
radial fall-off and, in practice, just like the result 
shown in Fig. ((SJi) for ^Pq in tetrad 7i, we are not able 
to obtain this radial behaviour numerically, and we end 
up getting just numerical error. 

In Fig. JSJ we show the comparison of with the 
news function; here the results are in better agreement 
but we still have no correspondence, the reason for this is 
to be found in the boost transformation, in fact a trans- 
formation like the one written in Eq. H21|l changes the 
value of ^'o according to 



(22) 




r (BH units) 



FIG. 7: (a) Top panel: the comparison of ^'4 (dashed line) 
calculated in tetrad Ti and — 7,t)u (solid line) for v = 80. (b) 
Bottom panel: the value for ^'4 alone. 




This leads us to the final conclusion that, in the linearized 



5 10 15 20 25 30 35 40 45 50 
r(BH units) 

FIG. 8: From top to bottom: the comparison of ^'o computed 
in the tetrad 72 and — 7,t)u for the time values of = 10, 
v\ = 40, D2 = 60 and vz — 80. The dashed line is 'i'o while 
the solid line is — j^w 

regime, the following relation must hold: 

We test this conclusion in Fig. where we have plot- 
ted the value of (^^^)^ (*o)t2- 
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FIG. 9: From top to bottom: the comparison of j4~^5'o 
(dashed line) calculated in tetrad T2 and —'y,vv (solid line) 
for the values of vo = 10, vi = 40, V2 = 60 and vs = 80. 



We emphasize that all the results that we have ob- 
tained in all the tetrads have wave-like profiles, although 
only one is the correct wave contribution. In practical 
numerical simulations one should really make sure that 
the tetrad in which the scalars are computed is a quasi- 
Kinnersley tetrad, otherwise the results, even if wave-like 
shaped, could be wrong. 



D. Energy calculation 

Having made sure that '^4 calculated in tetrad Ti is 
related, in the linear regime, to the Bondi news function, 
we can use its expression to calculate the energy radiated 
from the black hole. In section Hvl we have shown that 
the expression of the energy flux per unit solid angle is 
given by 



dvdft 



2^2 



47r 



(24) 



where we denote with $ the generic expression for the 
news function, being it 7_t, or J {'i'4)^j^j, dv. We can in- 
tegrate the expression in Eq. (|24|) on a 2-sphere in order 
to obtain the energy flux. For the sake of simplicity we 
take a sphere of radius ro , getting the result 



dv A-K J 




r (BH units) 

FIG. 10: From top to bottom: energy contribution of the 
I = 2 harmonic initial data. The three graphs represent the 
outgoing energy contribution for the values of i = 2, 4, 6. In 
each graph the two curves represent the value for the energy 
calculated using the Bondi news function (upper curve), while 
the lower curve uses the value of ^4 in tetrad Tl , as a function 
of the position of the observer. It is evident that at late times 
there is a convergence of the two values. This convergence 
seems to be less evident for the / = 6 case (lowest graph), but 
we expect this phenomenon to be purely numerical, because 
the numerical error on this multipole component is very high. 



The computation of the energy flux was the goal of 
|l6| . where the news function was used to calculate the 
amount of energy which is carried away by each spin- 
weighted spherical harmonics of the outgoing radiation. 
We can perform a similar calculation using ^'4 and com- 
pare our results, in order to have a further demonstration 
of the validity of our approach. Given the results de- 
scribed in Section lV Cl it is evident that also these results 
will be in good agreement, however we want to highlight 
their validity and to show their dependence on the posi- 
tion of the observer. 

Since we are interested in the energy contribution of 
each spin weighted spherical harmonic, we first have 
to perform the decomposition of the signal into spin 
weighted spherical harmonics contributions. This is done 
by introducing the quantity defined as 



$i (w, r) = 2TT J $ {v, r, y) 



(25) 



2io{y)dy, (26) 



where y = — cos and 1^2/0 is the spin weighted spherical 
harmonic of spin 2. 

Using Eqq. H25|l and (|25|l we can get an expression for 
the total energy emitted in each angular mode after the 
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FIG. 11: From top to bottom: energy contribution of the I = 3 
harmonic initial data. We show the two graphs corresponding 
to the dominant terms Z = 3, 5 in the emitted gravitational 
signal. Again here we compare the result coming from the 
Bondi news function with the one using 'I'4 in tetrad Ti . The 
results are similar to the ones shown in Fig. IIUII . 



evolution to a final time T, given by 



Ei{T)^^j^ [^i{v,r = r^)fdv. (27) 

We have performed some numerical simulations where 
the energy, using both the Bondi news function and 5*4, 
has been calculated. The results are shown in Figg. H10() 
and l(TT|l. 

Fig. HlOfl shows the result for a numerical simulation 
where the initial profile of the 7 has been chosen to be 
quadrupolar, i.e. using the spin-weighted spherical har- 
monic with I — 2, m = Q. The non-linearity of the prob- 
lem is translated into the fact that the evolution excites 
higher order multipolar terms. However, simmetry con- 
siderations allow only even multipolar terms to be ex- 
cited. In the picture we show the energy at time w = 80 
for the Z = 2,4, 6 terms. Such energy is calculated vary- 
ing the position of the observer and it is evident that, as 
soon as we push the observer further from the source, the 
two energy calculations coincide. On the other hand, nu- 
merical errors become stronger when going higher order 
multipole terms, which explains the not-perfect conver- 
gence for the I = 6 terms. 

Fig. Hll|) shows a similar simulation for an initial data 
with I — i, m = Q. Here again we expect the non- 
linearity to excite the other harmonics. Differently from 
the I — 2 case, we don't expect to have forbidden modes, 
however our numerical results show that the highest am- 
plitude modes are the odd ones, so we show only those 
modes. Anyway the results in this case are qualitatively 
equivalent to those obtained in the case of quadrupolar 
initial data. 



VI. CONCLUSIONS 

The problem of correctly extracting the gravitational 
signal in numerical simulations is of primary impor- 
tance. We believe that the use of the Weyl scalars 
of the Newman-Penrose formalism offers a promising 
method for wave extraction, as it applies to any formula- 
tion of Einstein's equations and, even more important, 
to any kind of background we end up with, being it 
Schwarzschild or Kerr. However the problem of iden- 
tifying the tetrad in which one is to compute the Weyl 
scalars still awaits a full solution. Recent work P, Q 
shows how to make the important first step of identify- 
ing an equivalence class of tetrads, the quasi-Kinnersley 
frame, of which the desired quasi-Kinnersley tetrad is 
a member. However, the problem of isolating the right 
tetrad out of this set is still under investigation. 

In the present work we have considered a non-trivial 
numerical scenario, namely the evolution of a non- 
linearly perturbed black hole using Bondi coordinates, 
in order to show the importance of the tetrad choice for 
the calculation of wave related quantities. This particu- 
lar scenario is well suited for a practical demonstration 
of the problems one would encounter if a careful choice of 
the tetrad for the Weyl scalar computation is not done. 
We have in fact shown that the computation of the Weyl 
scalars in an arbitrary tetrad, chosen by brute-force us- 
ing mathematical packages like GrTensor, would lead to 
wrong results for \E'4, which is the quantity that typically 
is supposed to contain the (outgoing) gravitational wave 
degrees of freedom. This fact is evident in our case, where 
we have compared directly the results for the Weyl scalar 
in two different tetrads 7i and T2, using the Bondi news 
function in determining that the Weyl scalar correspond- 
ing to the quasi-Kinnersley tetrad is the right one. In 
Section El we have shown that only the quasi-Kinnersley 
tetrad Ti gives results in agreement with the news func- 
tion. 

Finally, we emphasize the importance of singling out 
an appropriate quasi-Kinnersley tetrad from the quasi- 
Kinnersley frame 0,0- For instance, the tetrad T2 after 
exchange of the £ and n null basis vectors is related to 
the tetrad 71 by a boost. This example indicates that ev- 
ery tetrad in the quasi-Kinnersley frame will give results 
for ^'o and ^4 that will show no contribution from the 
background, so that the wave-like shape of the scalars 
could lead us to the wrong conclusion of having the right 
outgoing gravitational signal. As we have shown, this 
conclusion could well be far from reality. 
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